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We consider the problem of evolving nonlinear initial data in the close limit regime. Metric and 
curvature perturbations of nonrotating black holes are equivalent to first perturbative order, but 
Moncrief waveform in the former case and Weyl scalar i()4 in the later differ when nonlinearities are 
present. For exact Misner initial data (two equal mass black holes initially at rest), metric pertur- 
bations evolved via the Zerilli equation suffer of a premature break down (at proper separation of 
the holes L/M « 2.2) while the exact Weyl scalar ip4 evolved via the Teukolsky equation keeps a 
very good agreement with full numerical results up to L/M « 3.5. We argue that this inequivalent 
behavior holds for a wider class of conformally flat initial data than those studied here. We then dis- 
cuss the relevance of these results for second order perturbative computations and for perturbations 
to take over full numerical evolutions of Einstein equations. 
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I. INTRODUCTION 

There is a revival of the interest on perturbation the- 
ory of black holes since the work of Price and PuUin |^] 
who put forward the close limit approximation. This ap- 
proach considers the final merger stage of binary black 
holes as a single, perturbed black hole. They studied 
the Misner problem, two equal mass black holes initially 
at rest, and compared their results with full numerical 
evolution of Einstein equations for various initial separa- 
tions. The impressive agreement, even for not so small 
separations, triggered several researchers to test these 
ideas for initial data representing boosted towards each 
other, single spinning plus Brill waves, and orbiting black 
holes [|. 

Abrahams and Price found that if one does not 
linearize Misner initial data, but directly extract the £- 
multipoles from the exact 3-metric, the close limit ap- 
proximation breaks down for much smaller separations; 
in a regime where the linearized approach precisely agrees 
with full numerical calculations. In this paper we shed 
some light on the origin of this surprising result. 

Price and PuUin used the Zerilli equation in the time 
domain to evolve perturbations around a Schwarzschild 
black hole. There is an alternative approach to deal with 
perturbations of black holes (even with net rotation) due 
to Teukolsky |^ based on the Newman-Penrose formal- 
ism. The two approaches are related and equivalent when 
one deals with first order perturbations [^j^]. Here we 
want to explore how they behave under nonlinear com- 
ponents present in the initial data. This study is relevant 
to test the idea Q of using perturbation theory at the fi- 
nal merger stage of binary black holes taking over Cauchy 
data from full numerical simulations that started when 
the two black holes were fully detached (for instance near 



the ISCO). Our results will also serve as precise analytic 
benchmarks for this case. It is also useful in order to reli- 
ably use second order perturbation technics when initial 
data do not come in an analytic form, but as a numeri- 
cal table, thus making a complicated if not an impossible 
task to disentangling each perturbative order. 



II. INITIAL DATA 

Misner |^ found a solution to the conformally flat, 
time symmetric initial value problem representing two 
black holes at rest separated by a proper distance L 
parametrized by /ig (see below). For /io < 1-8 {L/M < 
3.3) a common event horizon encompasses the system, 
and for fiQ < 1.36 {L/M < 2.5) a common apparent hori- 
zon appears. 

The Misner three-geometry takes the form 
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Here we identified R, the conformal space radial coor- 
dinate, with the Schwarzschild isotropic coordinate. 



R = f=i (yr -\-Vr-2M 
The conformal factor is given by. 
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This metric represents an asymptotically flat three ge- 
ometry with total ADM mass M. 

The proper distance L between the throats can be writ- 
ten as 



(2.6) 



III. METRIC PERTURBATIONS APPROACH 

The theory of metric perturbations around a 
Schwarzschild black hole was originally derived by Regge 
and Wheeler for odd-parity perturbations and by Zerilli 
for even-parity ones. The spherically symmetric back- 
ground allows for a multipolc decomposition, labeled by 
(^m), even in the time domain. Moncricf Q has given 
a gauge-invariant formulation of the problem in terms of 
the three-geometry perturbations. For the head-on colli- 
sion of two black holes, that concern us here, only even 
parity modes are present. They are described by a single 
wave (Zerilli's) equation 
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Here r*^r + 2A/ln(r/2Af — 1), and the potential 
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where A=(f + 2){i — l)/2. The Moncrief wave function 
Q{irn)^ in terms of the metric perturbations in the Regge- 
Wheeler notation, is 
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Price and PuUin have consistently worked in the 
first perturbative order by linearizing Misner data and 
evolving them with the Zerilli equation. The excellent 
agreement with full numerical results up to values of /io ~ 
1.8 was somewhat unexpected. 



To linearize the 3-metric (2.1) we can use the pertur- 
bative notion of <i>^ — 1 <C 1. In this case, after expansion 
into Legendre polynomials, we get 



$4 



where 



M 
2R 



[M /R)'+^Pf (cos 0), 

=2,4... 

(3.4) 



OO 



(coth n^oY 
sinh nfiQ 



(3.5) 



Note that the expansion parameter, Mcoth(^o)/(4i?Si), 
has to be less than 1 for the above expansion to be 
valid. This condition is always fulfilled outside the ef- 
fective single hole event horizon, located at i? = M/2, 
for ^0 < 1.56865. For bigger values of no one has the 
inverse expansion analogous to that made in Eq. (2.21) 
ofRcf. 0. 

Since to first order we only have quadrupolar {£ — 
2) contributions, expansion ( |3.4[ ) has the nice feature of 
separating, in the common factor K2, all the dependence 
on the initial distance between holes. Thus allowing to 
compute the whole family of evolutions with only one 
integration of the Zerilli equation (3T), and then rescale 
waveforms by their corresponding K2(/io). In the Regge- 
Wheeler notation one finds that the o nly non-vanishing 
components of the perturbed 3-metric ( |2.l| ) are H2 — K. 
In this way one builds up the initial value of (with 
dtQ^ = 0) and then evolve to obtain waveforms, spectra, 
and radiated energies 

It was first found by Abrahams and Price |^ that if 
one does not linearizes the initial data, but extracts its 
multipoles from the exact metric as 
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and evolves using the Zerilli equation, one obtains for 
the radiated energy the disturbing results reproduced in 
Fig. |]. While for small separations, L/M < 2, corre- 
sponding to /io < 1.0, the agreement with full numerical 
computations |ll| and linearized perturbations is good; 
the energy curve presents a local maximum at fiQ « 1.25 
and then decreases to a local minimum at fiQ « 1.41. At 
this point the system radiates an order of magnitude less 
than the evolution of linearized initial data. For larger 
values of /xq the radiated energy rapidly increases and 
soon after /io ~ 1.49, when it crosses upwards the linear 
prediction, it overestimates the total energy radiated by 
several orders of magnitude. All this happens well before 
the linearized theory begins to deviate from the full nu- 
merical results {fiQ > 1.8). To unravel this paradoxical 
result we will first compute the corresponding radiated 
energies in an alternative formulation of the black hole 
perturbations. 
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FIG. 1. The total gravitational radiation computed in the 
close limit approximation via the Zerilli equation. The solid 
curve represents the evolution of exact Misner initial data 
for two equal masses black holes. At the bottom of the dip 
(L « 2.58) the system radiates one order of magnitude less 
than that predicted by the evolution of linearized initial data 
(dotted line). For comparison are given the results of full 
numerical integration of Einstein equations. 



IV. CURVATURE PERTURBATIONS 
APPROACH 

There is an independent formulation of the perturba- 
tion problem derived from the Newman-Penrose formal- 
ism 1^] that fully exploits the null structure of black holes 
allowing to uncouple for a single wave equation to de- 
scribe perturbations around Kerr black holes. The outgo- 
ing gravitational radiation is fully described in this gauge 
(and tetrad) invariant formalism by the Weyl scalar 



(4.1) 



where and ^fi^^ (together with and m^^) form the 
tetrad that span the spacetime. 

ipA fulfills the Teukolsky equation, which for the 
Schwarzschild case reads 
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coefficients cr = = A) and fix the form of m'' and its 
complex conjugate under transformations of type III 
(boosts) in such a way that the spin coefficient e = 0. 
Our tetrad is then 
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Using formulae (3.1) and (3.2) of Ref. ||T^ that give 
'04 and dt^'A in terms of the 3-geometry and the extrinsic 
curvature (here vanishing) we obtairj^ 
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The evolution of this initial data via the Teukolsky 
equation produce the results shown in Fig. ^ for the 
total energy emitted as gravitational waves. There is no 
dip at any value of the separation of the holes. The pre- 
dicted energy agrees with the full numerical results for all 
^0 < 1-8. For larger values of the separation (although 
we evolved here exact initial data for two distinct black 
holes) we overestimate the radiated energy at practically 
the same rate as the original Price-PuUin result who ex- 
trapolated the linearized, i = 2, initial data to values of 
(^0 > 1-58) beyond the radius of convergence of their 
expansion parameter, i.e. Mcoth(/io)/(4-REi) > 1. Cor- 
recting for this fact produce the £ = 2 piece of the curve 
labeled as "Linear" . 

How robust is this result? To answer this question we 
fist considered different tetrad rotations of type III to de- 
fine the relative normalization of the vectors P and 
while keeping fixed its null directions. The total radi- 
ated energy was quite insensitive to this choice. Some 
choices slightly improved the agreement with full numer- 
ical results, but we found not a priori justification for 
such choices. We thus keep tetrad (4.3), which in the 
Schwarzschild limit reproduces the background tetrad jj] 
used to write down the Teukolsky equation in Boyer- 
Lindquist coordinates. 

We also essayed two possible definitions of perturba- 
tive notion in the Newman-Penrose formalism. The first 
one, that we called perturbative, considers the exact Weyl 



where A = r'^ — 2Mr. 

The first step towards building up the initial -04 and 
dftpA is to find an instantaneous exact tetrad compatible 
with the data (^.1|). We have found it by choosing the 

and 71^ that generate shear free null congruences (spin 



*Note that Eqs. (3.1) and (3.2) of Ref. ^ are exact for 
our data if we just drop the (0) and (1) labels. Note also an 
obvious misprint there: The second addend in Eq (3.1) and 
the third one in Eq (3.2) should carry a factor 8 instead of 4. 
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FIG. 2. The total radiated energy as computed using the 
Newman-Penrose-Teukolsky approach. The solid line is the 
linear evolution of the exact initial data for two black holes 
initially at rest. There is complete agreement with full nu- 
merical results for L < 3.3 and with the Price-PuUin results 
(labeled as PP, £ — 2) all the way up to much bigger initial 
separations, keeping under control all the initial nonlineari- 
ties. Also shown two other notions of perturbative initial data 
in the Newman-Penrose formalism. 



tensor contracted with the background (Schwarzschild) 
tetrad 
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and a second possibility, that we caUed linear because 
it considers only Hnear terms in the conformal factor $, 
gives 
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The resulting energy is essentially unchanged in the 
fiQ < 1.8 regime and grows steeper than the exact ini- 
tial data for larger initial separations of the holes. This 
is because for the initial perturbative ip^s higher £ contri- 
butions are less under control than for the exact initial 
■04 ■ This is seen in the waveforms extracted far away 
from the system. While for the perturbative choices of 
■04 waveforms look fine up to /i « 1.6 those evolved from 
the exact ip^ reach fi « 2.0 without much higher £ con- 
tent. No dip was found at any value of the separation of 
the holes and in fact we could not reproduce Fig. ^ re- 
sults within the Newman-Penrose-Teukolsky formalism. 
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FIG. 3. The £ = 2-multipole projection of the exact Misner 
data is used to build up the Moncrief-Zerilli initial waveform. 
Waveforms are normalized by K2 in order to show only de- 
viations from the linear regime. Also shown is the Zerilli 
potential to stress that initial data only differ well inside the 
potential barrier. 



This gives us a measure of a certain robustness of this ap- 
proach against nonlinearities included in the initial data. 



V. A CLOSER LOOK AT THE PROBLEM 

Since both evolution equations, Zerilli's and Teukol- 
sky's are equivalent in the linear perturbations regime we 
are working in, the different results obtained must be re- 
lated to how -04 and Q handle the nonlinearities included 
in the Misner data. We then first plot in Fig. ^ the Mon- 
crief waveform normalized to ki such that in the linear 
regime all curves would superpose. Still for /uq = 1.0 the 
curve do not deviate much from that of the linear data. 
As we increase the initial separation of the holes we find 
an steady increase of the amplitude in the region r* < 
while for r* > the curves lie on top of each other. Most 
of the outgoing radiation is generated around the top of 
the Zerilh potential Vt at r*/2M « 0.95. The form of 
the initial data in terms of Qi do not show any particu- 
lar feature around /iq = 1.41. Their relative form contin- 
ues to grow in amplitude very close to the horizon, but 
remains almost unchanged further outside. We checked 
that the effect of the change of the expansion parameter 
at /iQ — 1-56865 described above is not responsible for 
the dip at /iq = 1.41. 

The quantity directly related to radiation dtQ initially 
vanishes for Misner data (time symmetric), in Fig. |^ we 
plot thus d\Q at i = 0. Observe again the good super- 
position for r* > 0. There is also a big depression for 
linear data with a minimum at r*/2M « ~3. As the 
two black holes start more separated and increases, 
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FIG. 4. The initial second time derivative of the Moncrief 
waveform normahzed by K2 to compare with the linear ini- 
tial data. While there is a good superposition of data form 
r* > 0, notable differences appear for r* < 0. the maximum 
amplitude decreases steadily for < < 1-41 and increases 
for /lo > 1.41. Also the location of the maximum ampli- 
tude monotonically recedes towards more negative r* as 
increases. This provokes a decrease of the frequency for which 
destructive interference of the outgoing radiation occurs. 

the amplitude of this depression decreases and its loca- 
tion recedes towards more negative 7'*'s. This continues 
so until we reach fiQ « 1.41. Then the amplitude begins 
to grow while still the location of the depression recedes 
towards negative r*'s. Is this relatively small decrease of 
the amplitude (around 20%) responsible for one order of 
magnitude less radiated energy? Actually since the de- 
pression is located very close to the horizon, well inside 
the potential barrier, most of the radiation generated is 
swallowed by the black hole and just a very small piece of 
it reaches infinity. To answer the above question one has 
to take into account the wave nature of the gravitational 
radiation. We know that a great deal of the radiation 
coming out to infinity is generated around the maximum 
of the ZeriUi potential at r*^a^/2M « 0.95. There will 
also be a piece of the disturbance generated at around 
the pick of d^Q. these two pulses will be out of phase 
by (2Ma;)(Ai/2M) where At is the time the pulse gener- 
ated close to the horizon takes to arrive at r* /2M k, 0.95. 
Depending on their relative phases these two pulses can 
produce constructive or destructive interference. Assum- 
ing At ~ Ar*, where Ar* is the distance between the 
pick of d^Q and the maximum of the ZeriUi potential. In 
the linear regime, this analysis show that the destructive 
interference appears for frequencies 2Ai"w > 1.2; too high 
to infiuence the total energy radiated at infinity. As we 
increase /io, the pick in Fig. ^ moves to the left, thus 
generating destructive interference at lower frequencies. 
In Fig. H we show that this effect is clearly visible in 
the three spectra for /io = 1.38, 1.41, 1.45. There is 
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FIG. 5. The spectrum of the gravitational radiation far 
away from the system as evolved by the ZeriUi equation. The 
solid line labeled with the separation parameter /io = 1 gives 
essentially the linear regime (we rescaled its amplitude by 
one fifth), we see how the suppression effect sets up sharply 
around /io = 1.41. For /to = 1.38, 1.41, 1.45 destructive in- 
terference occurs around frequencies 2Mu ~ 0.95, 0.83, 0.65 
respectively. 

destructive interference around 2Mw « 0.95, 0.83, 0.65 
respectively. The strong suppression of the radiation at 
/io = 1.41 is then due to destructive interference right at 
the frequency of the maximum of the linear spectrum, at 
2Mw ~ 0.8. 

By the same mechanism we can explain the sudden 
increase in the radiated energy for /io > 1.4. It is now 
the effect of destructive interference acting at too low 
frequencies and constructive interference at higher ones 
together with a dramatic increase of the amplitude of the 
initial data with respect to the linear regime. Finally for 
/io > 1.56865 the effect of the change in the perturbative 
parameter makes things blow up. 

Alternatively one can look at the initial form of '\\)^ and 
its time derivative cornputed via the linear transforma- 
tions given in Refs. in terms of the initial Moncrief 
waveform -0. This transformation can be seen as part of 
the process of linearization that takes place during the 
evolution with the ZeriUi equation. We observe then in 
Fig. U the same receding of the picks as the initial sepa- 
ration of the holes increase, and its separation from the 
maximum of the ZeriUi potential agrees very well with the 
prediction coming from the interference patterns seen in 
Fig. |. We wiU see that this behavior is in contrast with 
that of the direct computation of the exact ^"4 and 9* '04. 

VI. DISCUSSION 

Unraveled the mechanism that generates the Abra- 
hams -Price dip, there remains the question of why 
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FIG. 6. Initial and dtip4 derived from the Moncrief 
waveform 02 via Eq. (2.9) of Ref. § for different initial sep- 
arations of the holes. 
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FIG. 7. The exact initial i/)4 and dtip4 for the Misner geom- 
etry. The effect of increasing /xq is to decrease the amplitude 
and shift the maximum towards less negative r*. Curves tend 
to superpose in the region more relevant for radiation, r* > 0. 



nothing like this happens when we use the Newman- 
Penrose-Teukolsky approach to Schwarzschild perturba- 
tions. Again, a plot of the initial data for different sep- 
arations leads to the answer. In Fig. ^ we observe that 
the maximum of the initial data for increasing fio shifts 
toward increasing r*'s instead of more negative ones as 
happened for the Zerilli formalism. Thus regardless of 
a small decrease in the amplitude the destructive inter- 
ference occurs at too high frequencies to influence the 
outgoing radiation. 

The study of exact initial data evolved linearly has 
taught us several lessons: i) How the initial Misner data 
vary with respect to their linearized version and which 
aspects of them are relevant for computation of the grav- 
itational radiation that reaches infinity, ii) The inter- 
ference effects that occurs when one evolve these initial 
data. Note that the occurrence of interference is in gen- 
eral independent of nonlinearities. In fact it is rather sen- 
sitive to the shift of the maximum that they produce in 
the particular case of Misner data, iii) Since the Zerilli- 
Moncrief waveform and V'4 are two intrinsically different 
objects, they respond differently to nonlinearities. We 
have seen that for the initial data studied in this paper 
■04 is a better suited quantity to evolve linearly. 

There is still another way of evolving Misner data, 
namely nonlinearly, by direct full numerical integration 
of Einstein equations. In this way we of course evolve 
nonlinear initial data and obtain the correct behavior. 
In Ref. ijlj] we used the code Cactus to evolve the full set 
of Einstein equations for a single distorted black hole. We 
have found that even for small linear initial distortions 
we needed to solve the initial value problem nonlinearly, 
otherwise errors in satisfying the constraints would gen- 
erate numerically unstable evolutions. 

Our results are relevant to the idea of evolving fully 
numerically binary black holes starting from large sepa- 



rations (starting from Post-Newtonian initial data) un- 
til a common horizon encompasses the system and then 
let perturbation theory to take over The pertur- 
bation taking over full numerical method has the ad- 
vantage of optimizing supercomputer resources, concen- 
trating them in the region where the two black holes 
are completely detached and full nonlinear relativistic 
effects take place. Once a common horizon forms one 
can assume the close limit approximation to hold and 
continue the evolution with a single wave equation on 
the background of a Kerr black hole. It is very for- 
tunate that is the Newman-Penrose-Teukolsky approach 
rather than the Regge-Wheeler-Zerilli one that has this 
nice behavior in response to nonlinear Cauchy data since 
the Teukolsky equation can be generalized to rotating 
black hole backgrounds while the metric perturbation ap- 
proach do not. Unless there is a dramatic progress in full 
numerical technics, this perturbative-fuU numerical hy- 
brid, is our only chance to get waveform templates for 
inspiraling black holes before laser interferometers be- 
gin to operate. The effect described in Fig. [| inhibit 
us from using the Zerilli equation to that end. On the 
other hand, the Teukolsky evolution seems better suited 
for this marriage between full numerical and perturba- 
tive approaches. We also checked that the same general 
bad behavior of the Zerilli equation and healthy one for 
the Teukolsky equation holds when one considers Brill- 
Lindquist initial data (other solution of the conformally 
flat, time symmetric initial value problem). If this be- 
havior is also true in the astrophysically more appealing 
scenario of orbiting black holes is currently under inves- 
tigation. If we remain within the Bowen-York family of 
initial data (conformally flat and Longitudinal) we ex- 
pect a decomposition of the conformal factor of the type 

]l5[ <I> = ^Misner + ^reg, where ^reg IS prOpOrtioual tO 



6 



the square of the momentum of the holes and the square 
of their distance. This means that at least for black holes 
with small initial momentum the effects discussed in this 
paper should qualitatively still be present. 

Another situation where our results should be consid- 
ered is when one is interested in studying second order 
perturbations of rotating black holes . The perturba- 
tive approach to the Newman-Penrose formalism form a 
hierarchy of equations 



(N-l) 



(JV-l)i 



(6.1) 



where T stands for the background wave operator of 
the Teukolsky equation, V'l^'' is the waveform of the N 
perturbative order considered, and 5 is a source term 
formed by products of all perturbations of order lower 
than N. 

One can solve the above equations by evolving initial 
data order- by-order, successively reaching the next per- 
turbative stage. An alternative approach to that is to 
evolve exact initial data with the first order wave equa- 
tion (it has a vanishing source term in vacuum), and 
then evolve second order equations with vanishing initial 
data (expressing the source in terms of the first order 
perturbations). This is consistent to the desired second 
perturbative order. The justification to this approach 
can be found from a Laplace transform analysis. In the 
appendix of paper ||^ it is explicitly found that the initial 
data dependence can be included in the Teukolsky equa- 
tion as an additional source term, Sjz) and this holds at 
each order A'' in the above hierarchy of equations ( |6.1[ ) . 
Since the additional source term adds linearly and the 
wave operator is always the background one, summing 
over all A^ orders corresponds to evolve the first order 
equation with the exact initial data. Then one evolves 
all higher orders with vanishing initial data. 
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